Hadronic single inclusive k± distributions inside one jet beyond MLLA 



00 

o 
o 

c 

a 



43 ■ 

Or 
i ' 

Oh' 

<l> : 
43 . 



> 
On 
rn 

o 
o 



X 



Francois Arleo*, 1 Redamy Perez-Ramos, 2 and Bruno Machet 3 

'CERN, PH department, TH division, CH-1211 Geneva 23 
2 Max-Planck-Institut fur Physik, Werner-Heisenberg-Institut, Fohringer Ring 6, D-80805 Miinchen 
s Laboratoire de Physique Theorique et Hautes Energies , 
BP 126, 4 place Jussieu, F-75252 Paris Cedex 05 
(Dated: December 3rd 2007) 

The hadronic k ± -spectrum inside one jet is determined including corrections of relative magnitude 
O (y / a7) with respect to the Modified Leading Logarithmic Approximation (MLLA), at and beyond 
the limiting spectrum (assuming an infrared cut-off Qo = A qcd and Qo 7^ A QCD )- The agreement 
between our results and preliminary measurements by the CDF collaboration is impressive, much 
better than at MLLA, pointing out very small overall non-perturbative contributions. 
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Jet production - a collimated bunch of hadrons - in 
e + e~, e~p and hadronic collisions is an ideal playground 
for parton evolution in perturbative QCD (pQCD). One 
of the major successes of pQCD is the hump-backed 
shape of inclusive spectra, predicted in [1] within MLLA, 
and later discovered experimentally (see e.g. Q). Refin- 
ing the comparison of pQCD calculations with data taken 
at LEP, Tevatron and LHC will ultimately allow for a 
crucial test of the Local Parton Hadron Duality (LPHD) 
hypothesis [3[ and for a better understanding of color 
neutralization processes. In this Letter, a class of next-to- 
next-to-leading logarithmic (NMLLA) corrections to the 
single inclusive k x -distribution of hadrons inside one jet 
is determined. Unlike other NMLLA corrections, these 
terms better account for recoil effects and were shown to 
drastically affect multiplicities and particle correlations 
in jets [3| ■ We start by writing the MLLA evolution equa- 
tions for the fragmentation function D g (xl z, zEQ$, Qo) 
of a parton B (energy zE and transverse momentum 
k ± = zEQq) into a gluon (identified as a hadron h with 
energy xE according to LPHD) inside a jet of energy 
E. As a consequence of angular ordering in parton cas- 
cading, partonic distributions inside a quark and gluon 
jet, Q,G(z) — x/z Dq g (x/z,zEQo,Qo), obey the sys- 
tem of two coupled equations [5|] (the subscript y denotes 
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where &a( z ) denote the DGLAP [6| splitting functions, 
a 3 = 2ir/ AN c (3q(1 + y + A) is the one-loop coupling con- 
stant of QCD [13] and 

£ = (l/x) , y = \n(kjQo) , A = In (Qo/A qcd ) , 

(Qo being the collinear cut-off parameter), and where 

G = G(l)=xD%(x,Ee ,Qo), 
Q = Q{l)=xD%{x,EQ Q ,Q Q ). 

At small x <C z, the fragmentation functions behave as 

zEe Q ' 



B(z)^p h B (in J, In 



Qo 



-J =p h B (lnz + e,y), 



p B being a slowly varying function of two logarith- 
mic variables hi(zjx) and y that describes the "hump- 
backed" plateau |l| . In order to better account for recoil 
effects, the strategy followed in this Letter is to perform 
Taylor expansions (first advocated for in of the non- 
singular parts of the integrands in ()1I2|1 in powers of In z 
and ln(l — z), both considered small with respect to £ in 
the hard splitting region z ~ 1 — z = O (1) 

B(z) = B(l)+Bi(l)\nz + 0(\n 2 z) ; z^l-z. (3) 

Each ^-derivative giving an extra ^/al factor (see Q), 
the terms B t (l)\xiz and B e (1) In (1 - z) yield NMLLA 
corrections to the solutions of From and the 

expressions of the DGLAP splitting functions, one gets 
after some algebra (70 = 2N c o. s /-k) [8| 

Q(£, y) = 6(e) + ^ f d£' fdy'^' + y') (4) 
-£) + a 2 5(£' -£)M£\y')\G{£',y'), 
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G(£,y) = S(£)+ / d£' dy'^(£> + y>) 
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with ip e (i,y) = G e (£,y)/G(£,y). The MLLA coefficients 
5i = 3/4 and a% « 0.935 are computed in [B| while at 
NMLLA, we get Q: 
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Computing the NMLLA partonic distributions inside a 
quark and gluon jet, Q(z) and G(z), is the first step to 
determine the double differential spectrum d 2 N/dx d8 of 
a hadron produced with energy xE and at angle O with 
respect to the jet axis identified with the direction of the 
energy flow (see Q). As shown in it is given by 



d 2 iV 
dxdlnG 



(8) 



where F A is given by the convolution of two fragmenta- 
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f duDi o (u,E<d ,uEO)D h A (-,uE<d,Q 



(9) 

u being the energy fraction of the intermediate parton 
A. D Aq describes the probability to emit A with en- 
ergy uE off the parton Aq (which initiates the jet), tak- 
ing into account the evolution of the jet between Go 
and G. D A describes the probability to produce the 
hadron h off A with energy fraction x/u and trans- 
verse momentum k ± w uEQ > Qq (see Fig. [1]). As 




FIG. 1: Inclusive production of hadron h at angle Q inside a 
high energy jet of total opening angle Oq and energy E. 



discussed in [9(, the convolution © is dominated by 
u r-j 1 and therefore D Ag (u, EQ ,uE<d) is given by 
DGLAP evolution Q. On the contrary, the distribution 
D A = ID\ (z,uEe,Q ) = D\{l + \nu,y) at lowi:«« 
reduces to the hump-backed plateau, 



D A (l + lnu,y) X ^ U Pa (1 + \yiu,Y b +lnu), (10) 

with Yq = £ + y = ItlEQ/Qq. Performing the Taylor 
expansion of D to the second order in (In it) and plugging 



it into Eq. leads to 



xF Ao ~ E / duuDi o (u,Ee ,uEe)D h A (£,y) (11) 



j duu\nuD% a (u,E<d Q ,uEQ) 
- V f duuln 2 uDi o (u,EQ ,uEQ) 

2 A U 



dD h A (£,y) 



d£ 

d 2 D h A (£,y) 
d£ 2 



The first two terms in Eq. ([TT]) correspond to the MLLA 
distribution calculated in [9J when D A is evaluated at 
NLO and its derivative at LO. NMLLA corrections arise 
from their respective calculation at NNLO and NLO, 
and, mainly in practice, from the third line, which is 
new. Indeed, since x/u is small, the inclusive spectrum 
D A (£,y) is the solution of the next-to-MLLA evolution 
equations (|4]) and ([5]). However, because of the smallness 
of the coefficient a 2 (see (O), G(£,y) shows no significant 
difference from MLLA to NMLLA. As a consequence, we 
use the MLLA expression for G. It is determined here 
from a representation in terms of a single Mellin trans- 
form of confluent hypergeometric functions (see Eq. (24) 
of [n|), well suited for numerical studies [la]. The NM- 
LLA quark distribution Q{£, y) can then be deduced from 
G(£,y) using (g]) and ([5]), which yields 

Q(i,y) = ^f[G(l,y) + (ai-ai)Gt(l,y) (12) 



+ 



(ai (ai - aij + a 2 - a 2 ) Gu(£, y) + 0(jq). 



The functions and Fjf are related to the gluon distri- 
bution via the color currents (C) g . g defined as: 



xf: 



Nr. 



G(i,y). 



(13) 



(G)g.q can be seen as the average color charge carried by 
the parton A due to the DGLAP evolution from Aq to A. 
Introducing the first and second logarithmic derivatives 



of D'X, 



il)A,t{l,y) 



1 dD A (£,y) 
D h A (£,y) d£ 

1 d 2 D A {£ lV ) 
D A (£,y) d£ 2 



= 0(a s ), 



Eq. (jlip can now be written as 
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The scaling violation of the DGLAP fragmentation func- 
tion neglected in the last approximation is a 0(a s ) cor- 
rection to (u). It however never exceeds 5% [8j of the 
leading term and is thus neglected in the following. Us- 
ing (fig]), the MLLA and NMLLA contributions to the 
leading color current of the parton Aq = <?, q read 
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The MLLA correction, O (\/aTj , was determined in 0] 
and the NMLLA contribution, 0(a s ), to the average 
color current is new. The latter can be obtained from 
the Mellin moments of the DGLAP fragmentation func- 
tions 



&uui- x Di {u,£), 



leading to 



hi In 2 u)4„ 



—vi o (j,aEQ )-am) 
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Plugging (jTTJ) into fl6]). the NMLLA color currents for 
gluon and quark jets are determined analytically For 
illustrative purposes, the LO, MLLA, and NMLLA av- 
erage color current of a quark jet with Yq = 6.4 - cor- 
responding roughly to Tcvatron energies - is plotted in 
FigJ^] as a function of y, at fixed 1 = 2. As discussed 
in [9(, the MLLA corrections to the LO color current are 
found to be large and negative. As expected, the correc- 
tion O (a s ) from MLLA to NMLLA proves much smaller; 
it is negative (positive) at small (large) y. 

This calculation has also been extended beyond the 
limiting spectrum, A 7^ 0, to take into account hadroniza- 
tion effects in the production of "massive" hadrons, m = 



O(Q ) [10(. The NMLLA (normalized) corrections to 
the MLLA result are displayed in Fig. [3] for different val- 
ues A = 0,0.5,1. It clearly indicates that the larger A, 
the smaller the NMLLA corrections. In particular, they 
can be as large as 30% at the limiting spectrum (A = 0) 
but no more than 10% for A = 0.5. This is not surprising 
since A 7^ (Qo 7^ A cd ) reduces the parton emission in 
the infrared sector and, thus, higher-order corrections. 

The double differential spectrum d 2 N/dyd£, Eq. ([8|), 
can now be determined from the NMLLA color cur- 
rents |H]) using the MLLA quark and gluon distribu- 
tions Integrating it over I leads to the single inclusive 
y-distribution (or k ± -distribution) of hadrons inside a 
quark or a gluon jet: 



diV 
dy 



= k 



dN 
dk~ 



Yo -y 



d£ 



d 2 N 
didy 



9,1 
(18) 

The MLLA framework does not specify down to which 
values of I (up to which values of x) the double differen- 
tial spectrum d 2 N / dy d£ should be integrated over. Since 
d 2 N/dyd£ becomes negative (non-physical) at small val- 
ues of £ (see e.g. @), we chose the lower bound £ m i n so 
as to guarantee the positiveness of d 2 N/dyd£ over the 
whole £ m i n < £ < Yq range (in practice, £ 9 min ~ 1 and 
Anin ~ 2). 

Having successfully computed the single k ± -spectra in- 
cluding NMLLA corrections, we now compare the result 
with existing data. The CDF collaboration at the Teva- 
tron recently reported on preliminary measurements over 
a wide range of jet hardness, Q = EQq, in pp colli- 
sions at y/s = 1.96 TeV [12[. CDF data, including sys- 
tematic errors, are plotted in Fig. |4] together with the 
MLLA predictions of and the present NMLLA calcu- 
lations, both at the limiting spectrum (A = 0) and taking 
A = 250 MeV; the experimental distributions suffer- 

QCD r 

ing from large normalization errors, data and theory are 
normalized to the same bin, ln(fc_ L /lGeV) = —0.1. The 
agreement between the CDF results and the NMLLA dis- 
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FIG. 4: CDF preliminary results for the inclusive k ± distri- 
bution at various hardness Q in comparison to MLLA and 
NMLLA predictions at the limiting spectrum; the boxes are 
the systematic errors (their lower limits at large k ± are cut 
for the sake of clarity). 



tributions over the whole k ± -range is particularly good. 
In contrast, the MLLA predictions prove reliable in a 
much smaller k ± interval. At fixed jet hardness (and thus 
Yq ) , NMLLA calculations prove accordingly trustable in 
a much larger x interval. 

Despite this encouraging agreement with data, the 
present calculation still suffers from various theoretical 
uncertainties, discussed in detail in Q. Among them, the 
variation of A - giving NMLLA corrections - from the 

QCD ° ° 

default value A = 250 MeV to 150 MeV and 400 MeV 

QCD 

affects the normalized k ± -distributions by roughly 20% 
in the largest \n(k ± /l GeV) = 3 GeV-bin at Q = 100 
GeV. Also, cutting the integral (|18p at small values of I 
is somewhat arbitrary. However, we checked that chang- 
ing i 9 nlin from 1 to 1.5 modifies the NMLLA spectra at 
large k ± by ~ 20% only 16]. Finally, the k ± -distribution 
is determined with respect to the jet energy flow from 
2-particle correlations (which includes a summation over 
secondary hadrons), while experimentally the jet axis is 
determined exclusively from all particles inside the jet. 
The question of the matching of these two definitions at 
O (a s ) accuracy goes beyond the scope of this Letter. 

The NMLLA k ± -spectrum has also been calculated 
beyond the limiting spectrum, as illustrated in Fig. [5l 
However, the best description of CDF preliminary data 
is reached at the limiting spectrum, or at least for small 
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FIG. 5: CDF preliminary results (Q = 119 GeV) for inclusive 
k ± distribution compared with NMLLA predictions beyond 
the limiting spectrum. 



values of A < 0.5, which is not too surprising since these 
inclusive measurements mostly involve pions. Identifying 
produced hadrons would offer the interesting possibility 
to check a dependence of the shape of k ± -distributions on 
the hadron species, such as the one predicted in Fig. [5l 

To summarize, single inclusive k ± -spectra inside a jet 
are determined including higher-order O (a s ) (i.e. NM- 
LLA) corrections from the Taylor expansion of the MLLA 
evolution equations and beyond the limiting spectrum, 
A 7^ 0. The agreement between NMLLA predictions and 
CDF preliminary data in pp collisions at the Tevatron is 
very good, indicating very small overall non-perturbative 
corrections. The MLLA evolution equations for inclusive 
enough variables prove once more (see e.g. [(|) to include 
reliable information at a higher precision than the one at 
which they have been deduced. 
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